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We calculate azimuthal correlations in dAu collisions at different rapidities and centralities 
and argue that experimentally observed depletion of the back-to-back pick can be quantita- 
tively explained by gluon saturation in the Color Glass Condensate of the Gold nucleus. 



I. INTRODUCTION 



The subject of this paper is quantitative investigation of azimuthal two-particle correlations in 
deuteron (d) - Gold (Au) collisions at RHIC. These correlations are an important tool in studying 
properties of the Color Glass Condensate (CGC) [THS] - a coherent quasi-classical state of fast 
gluons and quarks discovered at RHIC [6]. In heavy ion collisions they are the main source of 
azimuthal anisotropy at semi- hard transverse momenta [^^S]. While the azimuthal correlations in 
heavy-ion collisions are smeared by the collective 'flow' effects, those in deuteron-heavy ion (dA) 
collisions allow for the direct observation of the CGC effect on the correlation function [lOj. A 
qualitative analysis of [TU] predicted most of the features of the azimuthal correlation function. 
However, a more precise data that recently has become available urges for a more quantitative 
approach, which we undertake in this article. 

The general form of the back-to-back correlations in hard processes in the kinematic region 
of large x follows from momentum conservation in the leading 2—7-2 process. In the center-of- 
mass frame the outgoing partons move in opposite directions with the same momenta implying 
a sharp pick at the azimuthal opening angle of A(p = vr independently of rapidity and centrality. 
It is essential for this argument that incoming partons are approximately on-mass-shell. This 
assumption breaks down when the produced particle transverse momenta become comparable with 
the characteristic transverse momentum scale of fast hadron or nucleus. At small x this scale is 
called the saturation scale Qs [IJ ^'^d is known to increase with centrality and decrease with x 
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[H [21 [TTl I12j . see (16). At the central rapidity at RHIC, its numerical value for a Gold nucleus is 
Qs — 1.5 GeV. Therefore, production of particles with comparable or smaller transverse momenta is 
strongly affected by CGC. Since x of nucleus decreases towards the deuteron fragmentation region, 
one expects even stronger effect of CGC on particle production in the forward direction |13H15] . 
Since the incoming partons are off-mass-shell, 2—7-3 process becomes the leading one in particle 
production yielding the uncorrelated background. This is natural, because CGC is a classical 
field in which correlations are suppressed. Correlations appear as quantum fluctuations around 
the classical field and correspond to 2 — )■ 4 processes. The corresponding azimuthal correlation 
function also picks at A(f) = vr, which, however, is not as sharp as in the hard processes since 
the momentum conservation no longer requires the outgoing partons to be back-to-back. This is 
the main source for depletion of back-to-back correlations for particles separated by small rapidity 
interval Ay <^ I. At large rapidity separations, the rapidity gap becomes filled with semi- hard 
partons further depleting the correlation. 

At present there are two heuristic approaches to take the CGC effects into account for observable 
quantities of interest. One is based on the 'dipole model' |16j in which one reduces the relevant 
scattering amplitudes to a product of light-cone "wave-functions" and combinations of CGC field 
correlators in the configuration space. These correlators satisfy a set of evolution equations with 
certain initial conditions. The advantage of this approach is that it rests on accurate theoretical 
treatment of the gluon saturation region. In this approach the double inclusive gluon [17J, quark- 
anti-quark |18f[20j and valence quark - gluon [2lJ cross sections were calculated assuming large 
rapidity gap Ay between the produced pair of partons (multi-Regge kinematics, MRK). It was 
argued in [8l [22] that this approximation misses important features of azimuthal angle dependence 
at not too large Ay. At the moment, extending the dipole model to the region of Ay ~ 1 presents 
a serious computational challenge. 

Another approach is 'A^j'-factorization', which assumes that 2 — )• n process and the two-point 
correlation functions of CGC fields can be factored out. In this approximation, the 2—7-4 ampli- 
tudes were calculated for an arbitrary Ay (quasi multi-Regge kinematics, QMRK) in p3ti26l [28] 
for gg — )• ggqq and in [221 E31 |27j for gg — )- gggg processes. The corresponding result is given by 



(10)-(13). Although generally fc^-factorization fails in the gluon saturation region, there are valid 
reasons to believe that it provides a reasonable approximation of the observed quantities. Indeed, it 
was proved that /c^— factorization provides the exact result for the cross section for single inclusive 
gluon production in the leading logarithmic approximation (LLA) Q |29] (though there is a sub- 
tlety in the definition of the unintegrated gluon distribution ip \15\ [29] ) . Although A:T-factorization 
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fails for the double-inclusive heavy quark production, the deviation from the exact results is not 
large at RHIC energies |30j . At transverse momenta of produced particles much larger than Qs, 
fcj'-factorization rapidly converges to the exact results. There are also numerous indications that 
/cT-factorization is phenomenologically reliable. Thus, the KLN model |31H34j that relies on kx- 
factorization provides successful description of experimental data. Recently, long range rapidity 
correlations in heavy- ion collisions were calculated in |35H38| . Although their theoretical results 
are more general, the phenomenological applications assume /cr-factorization. In view of these 
arguments, we will use /c^-factorization to compute the azimuthal correlations in this work. 

The paper is structured as follows. In Sec. [IT] we compute the azimuthal correlation function at 
small rapidity separations Ay ^ 1 at midrapidity y = and forward rapidity y = 3. In Sec. Ill we 
consider azimuthal correlations of particles separated by large rapidity gap Ay = 3. We discuss a 
possible effect of gluon evolution in the gap. Our calculations in both sections agree reasonably 
well with experimental data and provide an additional support for the CGC description of the 
nuclear wave function. We conclude in Sec. IIVI 

II. CORRELATIONS AT lyr - y^l < 1 

Azimuthal correlation function is defined as 

where dN/d{lS.(j)) is the number of pairs produced in the given opening angle and A'trig is the 
number of trigger particles. The number of pairs is given by 

dN ^ /" , /" , /" , /" , ( dNtrig (iiVass , dNcoTT \ 



27T [dkTkT fdyr fdkAkA ! dyAi "^^'"^ \ "^^'"^ ^ 

J J J J Xd^krdyT d^kAdyA d^hrdyr d'^kAdyAJ 



(2) 



d(A0) 

where k-p and yx are the transverse momentum and rapidity of the trigger particle and k^ and yA are 
the transverse momentum and rapidity of the associate one. We denote kx = etc. throughout 
this paper. The first term on the r.h.s. of ^ corresponds to gluon production in two different 
sub-collisions (i.e. at different impact parameters) and therefore gives a constant contribution to 
the correlation function, whereas the second term on the r.h.s. describes production of two particles 
in the same sub-collision. The number of the trigger particles is given by 

N,,;, = 2.jikrkrjiyr^^. (3) 
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Expression for the single inclusive gluon cross section is well-known |1H [2^ [59H45j . The corre- 
sponding multiplicity reads 

mUy = C7SI¥ J dqi^n{x+,q,)^Aix-,ik-^,) ), (4) 

where the unintegrated gluon distribution function ip is simply related to the gluon distribution 
function as 

xG{x,Q^)= / dq^^{x,q^). (5) 



x± are the fractions of momenta of incoming nucleons carried away by the t-channel gluons; in the 
center-of-mass frame 

x± = ^e^y. (6) 

Equation Q is derived in the multi-Regge kinematics (MRK) x± <^\. 
The correlated part of double-inclusive parton multiplicity is given by 

.2 /■ r ^2 



d^kTdvTd^kAdyA ~ ^^CTSI] ~W J ^ ^ 

X ipD{xi,ql) ^A{x2,ql) A{(\^,(1^,\t^^A^VT - Va) , (7) 



where 



kreVT + kAcy^ kre-^T + kAe'y^ 
Xi = , X2 = T= . (8) 

The amplitude A was computed in the quasi-multi-Regge-kinematics (QMRK) in [22l [271 EH] and 
recently re-derived in p3j (the gg — t- ggqq part was calculated before in |24H26j ) . In QMRK one 
assumes that xi,X2 <C 1, but Ay is finite. To write down the amplitude A it is convenient to 
introduce the kinematic invariants in a usual way. We will also denote kj^ = k^^, k^ = ^2 
Ut — VA = to simplify notations. We have [22] 

s = 2{kik2 cosh(Ay) — k^ • kg) , (9a) 

i = -kik2e^y , (9b) 

u = -{^^-^2? -kik2e-^y , (9c) 

S = X1X2S = kl + kl + 2kik2 cosh(Ay) . (9d) 



Now the amplitude A reads 



A — Agg^gg + ^-^^Agg^gq , (10) 
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where 



1 ^ 1 qhl 



e-^y 1 

+ TT0T0 + 



iu Aiuklkl 4ikik2 ^ukik2 Akfk"^ 



1 



2 / 11 \ 1 5] 

--{1 + kik2{-. - -) sinh(Ay) + ^r^{l + -) cosh(A?/)- 

S \ t U J 2KiK2 s 



1l 



4s ki t k2 



u 



.|i[(l + ^e-^.)l+(l + ^e^^)4] 
4s k2 t ki u 



1 1 ( ill - q,ni2 - - kfki 



+ 



til 



1 / fe - - kik2e-^y e\ / (ki - Si)' - ^1^26 



Ay 



4 Uk2-q^y + kik2e-^y s J Uk^-q^)'^ + kik2e 



Ay 



+ 



E 



with 



and 



E={%- %) ■ (ki - k2) - ^{ql - ql){kl - kj) + 2A;iA;2 smh(A2/) (^1 - ?L±M^ 



(11) 



(12) 



A 



gg^gg 



' sS 



1 + kxk2 sinh(Ay)(^ - ^) 
t u 



in J 



+ 



1 Aka - q^)' - fcifcae'^^ _ e\ / (k^ - q^)^ - ^1^26 ^^ ^ ^ 

2 Ife - q,)2 + feifeae-Af « ) I (k^ - q^)2 + fci^ae 



Ay 



+ 



+ 



[i^-cj^yn.-ci^y-ktk^ 

iu 



tu 



(13) 



In the Leading Logarithmic Approximation Ay ~ 1 /as ^ 1 the amphtude factorizes (Multi-Regge 
Kinematics) becoming 



— ■^gg^gg — 7.27,2 



Ay — )• 00 . 



(14) 



The double-inclusive cross section in MRK was first calculated in 

For numerical calculations we need a model for the unintegrated gluon distribution function (p. 
In spirit of the KLN model |3H - l33j we write 



27r2 



a,9 



)(1 



(15) 



At e"^/"" < X < 1 and q^ Ql (15) and ^ yield xG = InQ^ as required. The saturation 

scale is given by 



(16) 
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where xq = 3.04 • 10"^ and A = 0.288 are fixed by fits of the DIS data [Ml S?]. The couphng 
constant is fixed at = 0.3. 

It has been pointed out in (22] that due to 1 — )■ 2 gluon sphttings the double-inclusive cross 
section has a collinear singularity at s — )• 0, i.e. it is proportional to [(Ay)^ + (Acf))'^]^^. Such 
singularities are usually cured at the higher orders of the perturbation theory. Additional contri- 
butions to the small angle correlations arise from various soft processes including resonance decays, 
hadronization, HBT correlations etc. Because the small angle correlations are beyond the focus 
of the present paper we simply regulate it by imposing a cutoff on the minimal possible value of 
the invariant mass s. This is done by redefining the amplitude as A ^ As/{fP' + s). For each 
kinematic region, parameter jj, is fixed in such a way as to reproduce the value of the correlation 
function in pp collisions at zero opening angle Ac/) = 0. 

Expressions for the single and double inclusive cross sections written in this section are derived 
in the framework of /c^-factorization. As mentioned in Introduction, at the LLA this factorization 
gives the correct result for the single inclusive cross section. However, it fails for the double- 
inclusive cross section. The double inclusive gluon production was calculated in [17J and the double- 
inclusive quark- antiquark production in jl8H20j both in the LLA. Unfortunately, the LLA misses 
many important details of the azimuthal angle dependence. On the other hand, /cT^-factorization 
is known to give results that are in qualitative agreement with a more accurate approaches, but 
misses the overall normalization. Therefore, in order to correct the overall normalization of the 
cross sections we multiply the single inclusive cross section Q by a constant Ki and the double- 
inclusive one ([7|) by a different constant K2 [7]. The correlation function C depends on both 
Ki and K2. However, the difference Ca = C(A</)) — C(A(/)o) depends only on the ratio K2/K1. 
We choose in such a way that C(A0o) is the minimum of the correlation function. This is 
analogous to the experimental procedure of removing the pedestal |48j . The overall normalization 
of the correlation function K2/K1 — which is the only essential free parameter of our model - is 
fixed to reproduce the height of the correlation function in pp collisions. 

Integrations over transverse momenta of the trigger and associated particles can be simplified 
if we recall that because the transverse momentum spectra of both single-inclusive and double- 
inclusive cross sections are very steep, integration over the range of momenta is approximately 
proportional to the value of the integrand at the smallest value of momentum. Namely, if f{k) oc 
l//c" with n ^ 1, we have J^'"'"' dk f{k) ~ f{kmm) kmm / (^ ~ !)• We absorb the unknown constant 
(n— 1) in the overall normalization of C\. This is the approximation that we use in our calculations. 

The results of the numerical calculations are shown in Fig.[T|Fig.[2]and Fig. [3} In these figures we 
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observe suppression of the bak-to-back correlation in dAu as compared to the pp ones in agreement 
with the experimental data. In Fig. [3] we also see the depletion of the back-to-back correlation as 
a function of centrality. Note, that at the time of publication the precise centrality classes of the 
data shown in the lower row of Fig. [3] were not known. 




FIG. 1: Correlation function at the central rapidity. Kinematic region is 4 < pT < 6, 2 < pA < pr (all 
momenta are in GeV), i/t = 3.1, i/a ~ 3. Left (right) panel: minbias pp (dAu) collisions. Data from |48) . 
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FIG. 2: Correlation function at forward rapidities. Kinematic region is px > 4, 1.5 < pA < Pt (all momenta 
are in GeV), yr = 3.1, yA — 3. Left (right) panel: the minbias pp (dAu) collisions. Data from [49] . 



In addition to gg — ?• gggg and gg — t- ggqq processes that we took into account in this section, 
production of valence quark of deuteron gq^ — )• gq^gg gives a sizable contribution at forward 
rapidities due to not very small value of x associated with deuteron (x ~ 0.2 for pT = 2 GeV 
at y = 3). Contribution of this process to azimuthal correlations was analyzed in [21\ [SU] in the 
framework of the dipole model in MRK. However, the corresponding expression in A;7^-factorization 
in QMRK is presently unknown thus preventing us from taking it into account in our calculation.* 
In-spite of this we believe that the general structure of the correlation function as well as its 

* Authors of [171 121) discussed the process gq^ — >■ ggqv assuming the colhnear approximation to the unintegrated 
valence quark distribution |51| . 
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FIG. 3: Correlation function at forward rapidities. Kinematic region is px > 2, 1.5 < pA < Pt (a-U momenta 
are in GeV), yx = 3.1, i/a — 3. Upper left (right) panel: minbias pp (dAu) collisions. Lower left (right) 
panel: peripheral (central) dAu collisions. Note: centrality of the theoretical calculation may not coincide 
with the centrality of the data (the former is not yet known at the time of publication). Data from 




centrality dependence is not strongly affected by the valence quark contribution. We plan to 
address this problem elsewhere. 

To conclude this section, we would like to remark on a possible effect of fragmentation on 
the azimuthal correlations. Fragmentation is traditionally taken into account by convoluting the 
parton spectra ([4]), ([7]) with the fragmentation functions D{z). Particularly, hadron spectrum can 
be obtained from the gluon spectrum ([4]) as following \i2>\ I52j 

and similarly for the double-inclusive spectrum ([t]). Function F{k) encodes an effect of the DGLAP 



evolution at high k [131 152j. We are required to include this function because the integral in (17) 
extends to small values of z up to the kinematically allowed ^min- At small z's the high-A; tail of 
the parton spectrum, which is not accounted for in our model, plays an important role. Recall, 
that in this paper we use an approximation by which integrals over the transverse momenta are 
approximated by the value of the integrand at the lower limit of integration (see the discussion 



above in this section). Therefore, we can approximate the z-integral in (17) be the value of the 
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integrand at zq ~ ^(l), where zq corresponds to the maximal value of the integrand as a function of 
z. The net result of this approximation is a shift of the parton spectrum towards higher momenta 
k — )• k/zQ. We checked that our results shown in Fig. [l] - Fig. [s] are insensitive to such a shift 
in the region 1 < zq < 1/3. In order to set a limit on a maximal error that we make assuming 
this approximation, we performed an explicit integration over z as follows. Using the experimental 
data on the hadron spectrum [53J we set F(k) oc (1 — k/ko)~^, with n = 6 and ko = 1.2 GeV. 
For the transverse momenta of interest here, this effectively leads to an additional factor of z^ at 



small z in the integrand of (17). Fragmentation functions are taken from [54J. The corresponding 
modification of azimuthal distributions of Fig. [T] is exhibited in Fig. |4] We observe that even the 




FIG. 4: Effect of fragmentation on the azimuthal correlation function. Solid lines are the same as in Fig. [T] 
Dashed lines represent a conservative estimate of the fragmentation effect as discussed in the text. 

most conservative estimate leads to a minor modification of the correlation functions in agreement 
with the approximations that we made. Similar conclusion holds also for the correlation functions 
at forward rapidity. 

III. CORRELATIONS AT |yT - 2/a| > 1 



If the trigger and associated particles are well separated in rapidity so that \yT—yA\ ~ l/^s ^ Ij 



we can apply the MRK approximation (14) to the double-inclusive cross section. The result then 



factorizes into a product of two ladder rungs each given by the real part of the LO BFKL kernel. 
The corresponding formula is 



d?kT dyT (PkAdyA ■k'^CfS±_ k^k\ 



d qi (fD{xi,qi) fA{x2, (ky + k^ - q^) ) , 



(18) 



where xi^2 are defined in ([8j). The advantage of the MRK approximation is that it allows taking 
into account a possible multi-gluon production in the interval between yj- and y^- Unfortunately, 
in this approximation one also looses many features of the azimuthal angle dependence that are 
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important for description of the backward pick at intermediate Ay, see SecjlTj Thus, we are facing a 
dilemma: either to use formulas of Sec. Ill] that give precise dependence on Ay but neglect evolution 



in the gap, or to take into account the evolution as discussed below (see (19)-(23)) but in the MRK 
limit Ay — )• oo. At present, there is no approach that would interpolate between these limits at 
intermediate Ay relevant for RHIC. Therefore, in this section we will calculate the correlation 
function in two limits, compare our results with the data and try to learn which approximation is 
more phenomenologically important at Ay = 3. 



Eq. (18) does not take into account a possible gluon emission in the rapidity interval between 
yx and yA- This is important when \yi — y2\ > ^/ois and may be important for the experimentally 
measured forward-backward rapidity correlations. Evolution in between the rapidities of the pro- 
duced particles can be included using the the AGK cutting rules ^55j and the known properties of 
the BFKL equation |56j as 



d-kTdyrd^kAdyA = ^^2^^ ^ / ^ '^i / ^ ^M^i, ^i) V^^(:r2, (k^ + k^ - q,) ) 

xG{\q^-krp\,\q^-kA\,yT -yA) , (19) 

where G is the Green's function of the BFKL equation [56]. It can be written as 

oo 

G{qi,q2,y) = ^2cos(nqi • CI2) Gn{qi, q2,y) ■ (20) 

n=0 

Functions G„ are given by 

Gniqi,q2,y) = ^2 / dul^] e2^^^"M^ (21) 

where 

Xn{i^) = HI) - Re m{n + l)/2 + iu]} , (22) 



and ip{i') = T'{z)/T{z). One can evaluate the i/-integral in (21) in the saddle point approximation, 
which is valid when | In (71/52! ^ dsV- We have 



Gniquq2, y) - -(°) ^ exp | _ '""^^''^li, ] J , (23) 

2-K^qiq2 yasyi)2[[n + l)/2\] \j -2asy'il^2[{n + l)/2\ 

with '02 (i^) = r"(2;)/r(2;). The numerical values of poly-gamma functions for n = 1,2,3 are listed 
below 



Xo(0) = ln4, xi(0) = 0, X2(0) = -2 + ln4. (24a) 
V'2(l/2) = -14C(3), V2(1) = -2C(3), V'2(3/2) = -14C(3) + 16. (24b) 
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It is easily seen that Go has positive intercept, Gi has zero intercept and G2, G3, etc. have 
negative ones. Therefore, Go dominates at high y. However, n = term does not produce any 
azimuthal angle dependence. It is the sub-leading n = 1 term that is responsible for the azimuthal 
correlations at large ut — Va- Thus, we use an approximation G w 2[Go + Gi cos(qi • q2)] with Go 



and Gi given by (23). This is the same approximation that was employed in [10] 



Numerical calculations are presented in Fig. [s] together with the experimental data from |57|,I58]. 
We observe that the shape of the correlation function is better described by ([7])-([l3| in agreement 
with the observation of [3 [22] that finite Ay corrections to the MRK approximation are essential for 
description of azimuthal correlations. On the other hand, it seems that to explain the magnitude 
of depletion one also needs to include the small x evolution effects in the gap Ay = yr — VA- 
Obviously, a more accurate description requires additional theoretical investigation of the finite 
Ay corrections. Data from [57^ [58] also shows a significant isospin effect (not displayed here) that 
probably originates in the valence quark contribution not taken into account in the present work. 
This isospin effect obscures the CGC contribution and requires a detailed analyses that we plan 
to do elsewhere. At LHC one can get rid of the isospin effect by considering correlations at large 
rapidity gaps away from the fragmentation regions. 




FIG. 5: Forward-backward correlations. Kinematic region is pt > 2, 1 < p^i < pT (all momenta are in 
GeV), Ut = Va = 0. Left (right) panel: minbias pp (dAu) collisions. Solid lines: calculations with ([7j) 
(exact 2 — > 4 amplitude, no evolution between the trigger and associate particles). Dashed line: calculations 
with ( [l9| ) (MRK approximation of 2 4 amplitude, includes evolution between the trigger and associate 
particles). Data from [57l[58] (forward tt'' and midrapidity h^). 



IV. CONCLUSIONS 



We calculated the azimuthal correlation function in dAu collisions using the approach developed 
by us before 0[8]. The results are presented in Figs. T|2|3 andjs} We demonstrated that CGC is 
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responsible for depletion of the back-to-back correlations in dAu collisions as compared to those in 
pp ones at small rapidity separations - at midrapidity and forward rapidity - and at large rapidity 
separations. Our results quantitatively confirm earlier arguments of |1U| . 
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